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Abstract 

Starting from the Abrikosov-Ryzhkin formulation of the ID random potential prob- 
lem I find closed functional representations for various physical quantities. These func- 
tional integrals are calculated exactly without the use of any perturbative expansions. 
The expressions for the multipoint densities correlators are obtained. Then I evalu- 
ate the mean square dispersion of the size of localized wave functions. As a physical 
application of the method, I find the expectation value of the persistent current in 
mesoscopic ring with arbitrary magnetic flux <£. (For small 3> this problem has been 
solved by O.Dorokhov). The case when the random potential has finite correlation 
length is considered too. 
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1 Introduction and definition of the model. 



Anderson localization is acknowledged to be a fundamental macroscopic quantum phenom- 
ena. The localization manifests itself most evidently in one dimensional case. The essence 
of the effect consists in all the eigenfunction of the Hamiltonian 

n = ~ + u(x), (i) 

ax z 

to be localized wave packets providing the potential U(x) is a random function of x. (See the 
rigorous formulations, detailed discussion and bibliography in the book jL|.) This statement 
remains valid in the high energy limit considered in the present paper. 

The only quantities that can be calculated directly are various averages over an ensemble 
of potentials U(x). The measure of this averaging is reconstructed from the space correlation 
properties of a sample at our disposal. In the simplest case of the white noise statistics it 
takes the form: 

VUexpU^ ju 2 (x)dx^ , (2) 

< U(x)U(x') >= Dd(x-x'). 

Here (-L, L) is the interval, which our system takes up. P.Anderson has shown || 
that the difference from zero of the density-density correlator can be used as a criterion of 
localization of the state with energy E. The correlator can be expressed as 



p E (x,x') = lim <££(£ - E n )\m n {x)\ 2 \^ n {x')\ 2 >-- 
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= lim lim * < \G(x,x'\E + ie)\ 2 >, (3) 
where ^ n {x) are the eigenfunctions of 7i: 

n^ n (x) = E n m n {x) 

and G(x, x'\E + ie) is the resolvent of H: 

(n-E)G(x,x'\E + ie) = 6(x-x'), (4) 

Indeed, the continuous spectrum wave functions are of the order of 1/L 1 ^ 2 in every space 
point, and the sum over n gives effectively the factor L, so pe(x,x') ~ 1/L — > 0. For a 
homogeneous in the average potential the probability to find the state localized about a given 
point is ~ but the wave function ^ n (x) on its own does not depend asymptotically on 
L here. Hence only normalizable states contribute to pe{x,x') in the thermodynamic limit. 
(Assume boundary conditions in the endpoints of the interval (-L, L) provide the hermiticity 
ofH.) 

To calculate quantities like (f|) two approaches have been developed. The first one is so 
called "phase formalism." It allows one, in principle, to derive partial differential equations 
of the Fokker — Plank kind for various averages over the ensemble (0) H |J. (See for the 
review of advances M.) However, such an approach gives explicit results only in the regime: 

D 

Em « !> ( 5 ) 
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corresponding to the quasiclassical kinetics. Then another method can be used ||: extraction 
and summation of the infrared- singular terms of perturbation theory series (for review see 
H). The direct performance of this program requires sophisticated constructions and tedious 
computations. 

It has been noted in |7|, § that the sum of leading terms in the above mentioned per- 
turbation theory corresponds to some expectation values for spin 1/2 placed in a random 
magnetic field with Gaussian statistics. [] 

We present here derivation of this spin model somewhat modifying the line of arguments. 

For the wave function of the particle we assume the following boundary conditions: 

-fj-* n (x = -L) = * n (x = L) = 0. (6) 
dx 

The Green function (|J) can be expressed in terms of the solutions u(x),u(x) of the 
initial-value problems: 

(U - E) u = (ft - E) u = 0, (7) 

u'(-L) = 0, u{-L) = 1, u{L) = 1, u{L) = 0; 

„, m _ 1 J u(x)u(x'), x < x' 
' W I u(x')u(x), x' < X 

Here W is the Wronskian of the functions u and u: 

W = —u'(x)u(x) + u(x)u(x). (8) 

All the physical quantities of interest can be defined through one solution only, say, u(x) 
(see below). One can introduce for the function u(x) the "plane- wave components" v\(x) 
and V2(x): 

Vi(x) = e~ ikx (u'(x) + iku(x)) 
v 2 ( x ) = - e ikx (u'(x) - iku{x)) ,E = k 2 (9) 



so that Vi = (f2 = 0) for the plane wave propagating from right to left (from left to right). 
The equation (0) is equivalent to the following first-order matrix equation: 



(10) 



d_ / vi(x) \ _ U(x)/2ik, U(x)e~ 2ikx /2ik \ / v x (x) \ 
dx { v 2 (x) ) ~ [ -U(x)e 2ikx /2ik, -U{x)/2ik ) \ v 2 {x) ) 

and reduction 

-v 1 (x) = v* 2 (x). (11) 

It is seen from ( |10|) that the derivatives v\ and v 2 with respect to x are small along with the 
potential U(x). That is, vi(x) and v 2 {x) are changed slowly compared to exp(±ikx). Let us 
rewrite flTOf) in more compact notations: 

b = (iip(x)s z + ( + (x)s~ + C"(x)s + ) v. (12) 

Here 

v = ( V f\ ) 

V M x ) J 



1 Conceptually the same method was used in the paper 
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tp(x) = -U(x)/k, (^(x) = ±iU(x) exp(±2ikx)/2k, (13) 

s z = a z /2,s ± = (a x ± ia y )/2 are the usual spin operators and the dot denotes here and 
below the x-derivative. The formal solution of (O) can be written in the form: 



v(x) = T(x, —L)v(—L) 
I x 

T(x,-L) =Texp J (ip(t)s* + ( + (t)s- + (-(t)s+)dt\, (I I) 

where the sign implies the product is ordered along the interval (— L, L). 

Let us consider the expectation value of some functional of t>i(x),t>2(x). Expanding the 
-exponential (|T4]) and this functional in a series in the fields (^(t),C ± (t) and performing the 
averaging over T>U (x) we obtain the result as a series in integrals: 

dtdt' < (p(t)ip(t f ) >, f dtdt' < ( + (t)C(t') >, 



J dtdt' < (pty^lf) >, J dtdt' < C + (t)( + (t') >, J dtdt' < C(t)C(t') >, (15) 

over some domains of the order of L. 

In the last three expressions we integrate fast-oscillating functions. Therefore these inte- 
grals remain restricted in their values with increasing integration intervals and fall with an 
increase of energy. On their turn, integrals of the first two kinds correspond to the infrared- 
singular contributions and grow linearly with L. Thus to leave in the perturbation theory 
series the terms dominating in the large L limit the correlators: < C + C + >> < C~C~ >> an d 
< <f, ^ > should be neglected. It is equivalent to the assumption that the fields cp and C ± 
are statistically independent and the weight of T*^ - averaging is phase invariant. For the 
white noise statistics the corresponding integration measure has the form: 

D^PC^expj-^y (aip 2 (x) + ( + (x)C(x)) dzj , (16) 

where 

a= 2P' fl= 8- (17) 
We shall consider below the parameter a as an arbitrary ones. (It does not enter the final 
results.) 



The formulae (II) and (|1~6| ) were first proposed for the one dimensional random potential 
problem in the work [[?[]. Our presentation of it does not refer to the existence of the Fermi 
level. It allows one to suppose that Abrikosov-Ryzhkin model has some universal features 
relevant to the infrared behaviour. The model can be easily generalized to random potentials 
with finite correlation length. It may be usable to study spectral properties of operators 
which are not random in the strict sense (see Conclusion). 

The terms neglected in deriving of QHD, (Pp are smaller by a factor ~ l/(kL) compared 
to the ones kept. Thus this model can be applied to the study of mesoscopic systems (see 
section 5), since the inequality l/(kL) <C 1 for sufficiently large k is compatible with I > L, 
where / is the mean free path. 

The authors of the paper have used the formulae (|14|), fll6|) to obtain the conductivity 
of a one dimensional metal. Unfortunately, the calculations have being carried out there 



4 



by perturbation theory method lead to cumbersome constructions, which are inadequate 
to simple model. In the present paper I solve this Abrikosov-Ryzhkin model exactly with 
the help of functional integration method. On deriving the path integral representation I 
find the multipoint correlators of arbitrary powers of the density. With the use of these 
expressions I evaluate the mean-square dispersion of the size of localized wave function. As 
a physical application of the method I calculate the mean absolute value of the persistent 
current in a mesoscopic ring with an arbitrary magnetic flux $. (For small $ it has been 
found recently in ||.) In Conclusion I analyze the localization length dependence on the 
correlation length of the random potential. I discuss also a quantity that could play the role 
of the order parameter describing localization. 

2 Functional representation for averaged functionals of 

v(x) 

It is impossible to express v(x) as a functional of the fields <f(x), (^{x) explicitly. The same 
problem arises when one undertakes an attempt to write out a closed functional represen- 
tation for the partition function of quantum Heisenberg ferromagnet. It has been solved in 
the works [P^|-[P^[ and here we take advantage of the method proposed there. 



The ordered exponential T(x, —L) is defined by the equation 

t = {i V {t)s z + (+(t)s- + C(t)s+)T (18) 

and the initial condition: 

T(x = —L, —L) = 1. (19) 
Let us consider the operator given as a product of usual matrix exponential: 

/ X \ / X I t \\ 

T(x, —L) = exp (s + ip~(x)^ exp is z J pdt\ exp s~ J dtijj + (t) exp \i J pdr x (20) 

V -L j \ L \-L J J 

x exp (— s + ip~(— L)^j . 
Here ^ (x) , p(x) are some new fields. It obeys the equation: 

f = + 2i) + ^-)s z + ip + s~ + ~ ~ ^ + 0T) 2 ) S+ } (21) 
and the last factor in ( pO"D gives equality: 

f(-L, —L) = 1. (22) 
Thus, the change of variables in the functional integral over the measure ([16|) : 

iip = ip + 2ip + ip~ , 
{- = ip--ipip--iP + (ip-) 2 (23) 

brings the ordered exponential T(x, —L) to the form (p0|): 

T(x,-L) = f{x,-L), (24) 
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and allows us to obtain an explicit functional integral representation for any physical quantity 
to be averaged. (Parametrization of SL(2, C)-valued functions on two-dimensional space 
analogous to the (|23"D has been used also in the paper fl3[| ). To accomplish the change of 
variables in functional integral, much like the usual integrals, we need know the map ( p3|) in 
one direction only: from (</>, to (p, ijj^). The Jacobian J~[p, if)^]: 

V(pV( + VC = Jlp^^VpViP+Vij- (25) 



depends on the regularization of the map fl23|) and on the kind of condition imposed on the 
field The latter is necessary since there is first-order derivative of tp~ on the right-hand 
side of (p3|). The periodic boundary condition renders the map (^) irreversible. Following 
papers [fLlp , fll2l we consider the field if)~(x) as obeying an initial condition: 



^-(-L)=^o, (26) 



but, unlike JTTJ] , [f[2"| , the concrete value of ?/>o will be picked as the situation requires. 

The regularization of the map ([23]) is determined by the physical meaning of the model: 
the white-noise correlator (^) is to be considered as the limit of a smooth symmetrical 
correlation function. Any such a regularization of the 5-function gives for the correlators: 

t t 
<C+(t) / C"(0 dt>) = (C"(t) / C + (t') dt>) (27) 
o o 

the limiting value equal to ™ what corresponds to the extension of a definition of the step 
function 9{x): 

6(0) = 1/2. (28) 

The discrete version of the change of variables ( p3| ) providing the equalities fl27|) has the 
form: = ^(Q^n = p(t n ), . . . , n = 1, . . . , M, t n = -L + , h = § -> 0, M -> oo) , 

C = ^(C - V'n-l) - 2 iPn ^ + V»n-l) - (V>n + ^n-l) 2 > (29) 

C = C 

All the over-diagonal elements of the differential matrix of the map (|29|) equal zero. Then 
the Jacobian being the determinant of this matrix, is equal to the product of the diagonal 
elements only: 

/ . L \ 

J = const exp — / pdt . (30) 

V L J 

Making the substitution (|23| ) into the measure ( |i6|) with the use of the expressions ([25]), ( |30| ) 
we obtain the weight of averaging over the fields (p, ip ): 



NVpVip+Vtfj- exp (S(p, tj^Y 

L . L 

S(p^ ± ) = - f dx (ap 2 + - {1 + Aa)ipip + ^- - (1 + 4a)(^+^) 2 ) + - f dxp. (31) 



— L — L 
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Here N is a normalization constant depending on ah. 

In calculating of the Jacobian (|25| ) we were considering (p, ip^ and (0, as sets of 
independent complex variables or, in other words, as distinct coordinate systems in the 
whole space C 3M of the fields' configurations. The conditions 

Im V = 0,( + = ({-)*, (32) 

being from the outset embedded in the model specify the surface E in C 3M along which 
the differential form U(p f\U( + or T>p/\T>%p + /\T>%p~ is integrated. From the point 

of view of the coordinates set (p, ip ± ) the equation (|32"D for £ is implicit. According to the 
Cauchu-Poincare theorem the integration surface can be deformed in an arbitrary way in the 
convergence domain while an analytical function is integrated. There exists a continuous 
family of surfaces (homotopy) situated as a unit in "perturbative" convergence domain, 
which includes both the surfaces £ and the "standard" one E': 

E' = {/mp = 0,^ + = (ip~)*}. (33) 

The word "perturbative" means here that we check the convergence in every order of the 



perturbation theory expansion. (This homotopy is presented explicitly in the paper [ 10 ] . ) 
Thus, treating functional integral as the sum of perturbation theory series |TJ| we can replace 
the surface of integration £ by the standard one £'. 

However, to pass from E to £' the expressions being averaged (and not just the action) 
should be written in the form allowing the direct analytical continuation from the surface £. 
It means constructively that the definition of any physical quantity in terms of the matrix 
elements of T(x, —L) must contain no complex conjugations. 



3 The density-density correlator expression in terms 
of the functions v(x) 

The formula ([3D defines the correlator pe{x, x') in terms of the singular at e — > +0 part of the 
Green function G(x, x'\E + ie). When we use the representation (|8|) the singularity appears 
owing to zeros of the Wronskian W(E) on the real axis. Neglecting e in the numerator (|8|) 
and substituting 

W(E±ie) = W(E)±ieW'(E) (34) 
into the denominator we obtain: 

/ u 2 (x)u 2 (x') 



Pe{x,x 



] -S(W)\ ,x' > x. (35) 



\ \W'{E)\ v 7 

W(E) does not depend on x and, thus, we may put in (§) x — L: 

W = u(L) (36) 

Being in the product with 5(W) = 5 (u(L)) the solution u(x) is proportional to u(x). The 
proportionality coefficient is determined by the conditions (^). So, for an arbitrary functional 
T [u(x)} the following equality takes place: 



F[u(x)]5{u{L)) — T 



u[x) 



8{u{L)). (37) 
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According to (36) we can express W'(E) in terms of the derivative of u(x) with respect to 
the energy E: 

W'(E) = *^>. (38) 
The function g(x) = du(x)/dE obeys the equation 

/ d 2 ^ 



dx 2 



+ U(x) - E) g{x) = u(x), (39) 



and the initial condition: 

g{ x = -L) = g\x = -L) = 0. (40) 

The substitution g(x) = q(x)u(x) leads to the first-order equation for q'(x); its solution gives 
us: 

X dy V 

g(x) = u(x) J — ^- J dy^iyt), (41) 

—L U —L 

and: 

m^T\ s « L » = TTn\ 5 = J"'^ 5 ■ ( 42 ) 
\W'(E)\ \g(L)\ ) u \y)dy 

-L 

Thus the correlator Pe{%, x') can be written via u(x) as follows: 

PeM = ( u2 WfW 5(u(L))),x> > x. (43) 
\\u'(L)\ J u 2 {y)dy ' 

-L 

In the high-energy limit fl5|) we can get rid of the 5-function and obtain a simple formula 
for pe{x,x') in terms of slowly varying amplitudes Vi^(x). Indeed, in neighbourhood of any 
given point Xq the function u(x) can be written in the form: 

u(x) = u s i(x) sin(fca; + 5) (44) 

, where the envelope u s i(x) and the phase 5 vary only slightly over distances of the order 
~ 1/k. Let us average the expression (|43| ) over the interval AL of the right endpoints' 
positions of our "space" (—L,L): 

L+AL 



p E (x,x') = — J p E {x,x')dL (45) 

L 

- < AL < - = /. (46) 
k a 

(Here we introduce the standard notation I for the localization length.) In the thermody- 
namic limit the functions Pe{x,x') and pe{x,x') coincide. On the other hand, the value of 
u(x) in a given point, by the construction, does not depend on the right endpoints position. 
The integral in the denominator of is determined by the envelope u s i{x) only. The 
variation of L from L to L + AL does affect it asymptotically. The averaging (|45|) is suffi- 
cient only for the factor S (u(L)) /\u'(L)\. Since the conditions ( pf ) mean that u s i(x) can be 
considered as a constant in the averaging interval, we obtain: 

1 L+AL 1 1 

az I dL wm S{u(L))= ^m)' (47) 



s 



We can derive similarly the relationship between u 2 sl {x) and u 2 (x), in particular: 

2 L+AL 

<i( L )*XL I dLu2{L) - (48) 



Substituting into (f43|)-(f48D the expression of u(x) via v(x), neglecting the contributions 
vanishing in the k — > oo limit and keeping in the numerator of (|^) the "resonance" terms 
only, we obtain: 



p E (x, x') « p E (x, I'jw-j U ; u I ) ) jX > >x . (49) 



ui(L)v 2 (L) / v l {y)v 2 {y)dy l 

-L 

The "non-resonance" terms containing the oscillating factors exp ±2ik(x — x') will result in 
exponentially small in aL contributions and, thus, can be neglected. 



4 Functional integration for correlators of the density 
-density type 

The form of the expression (|49]) allows direct analytical continuation from the surface 
Vi = (^2)* over the functions Vi^x). It is not restrictive to put exp(ikL) = 1 in the thermo- 
dynamic limit. (Only if the formula (f£9|) has been obtained already!) Initial condition for 
v(x) takes the form: 

v(-L) = ik ( j J . (50) 



To find v(x) we substitute into (|HD the expression (|20|) for the evolution operator T(x, —L) 
picking the quantity ipo to be equal to 1: 

Vo=-0"(-^) = l- (51) 

It yields the equality: 

v(x) = exp 



/ 



y>*}n x) h ^ 



and the expression for pe(x,x'): 



,ip (x)ip (x') exp { — i J pdt — i J pdt + i f pdt\, 

Pe{x,x>) « / — ^-^ s — — J -\,x' > x. (53) 

ip~(L) f^ij-(y)expl-i J^pdtj dy 



Here the averaging over T> pD^Vip' is carried out with the weight (pl|) . To calculate this 
functional integral we employ a trick similar to the so-called "bosonization" in the field 
theory models [15]. Using the identity: 

exp (-^(p,^)) = / Xfyexp (-£(77, p,^)) , 
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2 L ' L 



^(ry,p,'0 ± ) = - J dx(^{l + 4a)r] 2 + ap 2 + tp + if- + {l + Aa)(2r]-ip)^ + tp-) + % - J dxp, 

-L -L 

(54) 

and the gauge transformation: 

^±(x) = x ± exp I ±(1 + 4a) / dt (2n - ip) I , (55) 



— L 

we get rid of the non-linear terms in the action. The Jacobian of the rotation (|55| ) is equal 
to 

j(2 V -ip)dt\, (56) 

V ~L J 



Jr = const exp 



where the regularization (29) is taken into account. The fields n and p enter the equation 
(|53|) via the combination 



J (2(1 + Aa)n - Map) dt 



only. It is natural to consider it as a new integration variable: 

£ = 2(1 + 4a)?7 - Map, 

£{-L) = (57) 
T>pT>n = const T>pT>^ 

Then the Gaussian Dp-integration can be done easily and we obtain the expressions for the 
measure 

L „ L 



fir. 2 f • £(L) I 
constViV X + V X -e W \-—j dx£ 2 J dx X + X~ ^~ (58) 



and for the quantity to be averaged: 



/ /\ 1 / X' (X) X - (g) exp (-fls) - gCgO + ^)) \ - /rnN 

p E (a;, x ) = — — ( 1 ) , x > x, (59) 



-L 



(The asymptotic equality in the limit (|5|) is assumed.) The initial condition for the field 
X~{ x ) follows from (|51~1): 

X'(-L) = 1. (60) 
It means that x~( x ) contains both the fluctuating part xj( x ) an d the regular one: 

X -(x) = l + x}(x), X ](-L) = 0. (61) 

The component x]( x ) does not contribute to Pe{x, x') because the conjugated field does not 
appear in the broken brackets in floDp. Thus, the only averaging over the field £(x) remains. 
Its weight has the form 



exp (-^) iV'^exp (-i- / dx? - ^ ] (62) 

V a -L 
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Here the normalization constant N' is determined by the quadratic in £ term of the 
action: 

iV'P£exp^-^| dx^ =1. (63) 

The factor exp (— aL/4) provides the equality < 1 >= 1 for the averaging over the entire 
measure (|62|). Thus we arrive at the following path integral for the correlator pe(x,x'): 



Pe(x,x>) = ^N'exp (-if) X_ o P£exp / ^ 2 - ^ 



x 



-l 



x exp - + { J exp (-£(*)) dt 

N' 



Arrka 



exp f-^) /dA / ^exp (-± J dt (£ 2 + Xe< + e~^~^ = 

= N' (Anka)- 1 x (64) 



x I I doda oxp ( } | VL exp I --f- | (// |£ 2 + e"^ - — I e - ^ -5 ^' 5 . 



+oo / / \ / L / • 

// dada' exp (^) / D£ exp / (V 2 

-do v ^(-L)=o-'^(L)=o- \ -L 



4 y 



The last equality has been attained by changing of variables: 

\ = e- a \^i-a' (65) 

and by separating the integrals over the values of in the endpoints t = L and t = —L. 
The final path integral in (B|) is of the Feynmann-Kac type [HJ and it is equal to the 
following matrix element: 

Pe(x,x') = exp ( — —\ (Anka)" 1 x (66) 



4 

e^ 2 \ exp (—(L - a/)#) e"«exp (-(x' - x)E) e~* exp (-(x + L)#) |e«/ 2 ) 



x 

with the Hamiltonian: 

a o9 1 

-fl? H 

2 € 2a 

The function e^ 2 increase when £ — > oo and, consequently, it cannot be represented as a 
linear combination of the eigenf unctions of ff: 



# = -"df + ^e- ? . (67) 



#/„(£) = V/,(0, /,(0 = -Vi/ sinh2vrz/K 2w (-e^ 2 ) , (68) 

2 7T \OL J 

Still, explicit solution of the corresponding evolution equation leads to the following asymp- 
totic relation: 

exp(-T^) e^ 2 — exp (^) T (£) = £ exp (^) ^ (V«/ 2 ) . (69) 
Here K^z) is the standard notation for the modified Bessel function. The function 

To(Oe" ? = -K X (-e~^) e< 
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in its turn can be expanded in terms of the complete set (|68[) . Thus, the correlator Pe{x, x') 
is equal to: 

p E {x,x') = ^^ ex P ^~ a ^g — j (T (O e ~ f l ex P (-\x-x'\H) |T (Oe _? ) = 

ex ( oAx — 3/ 1 \ 00 / 2 \ f 00 1 ^ 

exp I J J dvv sinh27rz/exp \-^-\x - x'\j < j dyyK x {y)K 2iv {y) > = 



ir 3 k 



- f - ("^) f ^ 1)' - - D <™> 

The formula (|70 ) up to the redefinition a/2 = Z _1 coincides with the well known result 
|I[,P|. This method allows us to compute the high-order correlators as well. For example: 
(xi < x 2 < ■ ■ ■ < x 2m , m > 1) 



2irkp < §' m \xi,x 2 , ■ ■ .x 2m -i,x 2m ) = 
= 2ixk lim (ES(E - E^^x^^^x,)^ . . . n {x 2m ^)\ 2 ^ n {x 2m )\ 2q ) = 

L^oo n 

2 2qm (qm — 1)! e 2qm ~ 1 k 

lim lim — — (\G(x 1 ,x 2 \E + ie)..G(x 2m _i } x 2m \E + ie)\ 2q ) 

l^oo e^+o {2qm — 2)1 



2nk ( 1 nu 2 %x 2j ^)u 2q {x 2j )\ {W'{E)) 



- 2qm+1 5(W(E))j 



(71) 



— 2gm+l 



u'(L))- 1 *(«(£)) 



2m 



( L 



-2qm+l 



EI (v 1 (x j )v 2 (x j )) q fv 1 {y)v 2 (y)dy) 



u'=i 



i\T' 



exp — 



( Vl (L)v 2 (L))-^ 

ah" 
T 



(2a) 2qm - 1 (2qm-2)\ 

x If dada'e^')! 2 J V^exp (-£ / <ft (£ 2 + e"*)) exp f- g g£( % )) 

a(x 2m - xi)\ 



exp 



(2a) 2 « r "- 1 (2gm- 2)! 
x W Q(£)e~ q ^\e~^ X2 ~ X1 ^ e~ q ^e~( xz ~ X2 ^ e~ q ^ . . . e~ q ^e~ ( - X2m ~ X2m ~ 1 ^\To(^)e~ q ^) 



71 



m+1 



r(g) 



2m-2 



T(2mg - 1) (r(2g)) 2 Vr(?+l/2), 



exp 



X 



oo oo 2rn> 1 / 2 \ 2?n 2 

x /. / n ^ exp ( -^Ax,- ) Vj 81^2^ pW(v 1 )pW( I / 2m _ 1 ) n ^ s +i; 



o o i=i 

where Ax,- and the functions P^ q \v) and Q^iy, i/) are defined as follows: 



1 



COsh 7Tf 



9-1 



2 , ^2 



(g-l/2) 2 + z, 2 jn((j-l/2) 2 + 

3=1 



Q iq) '(«/,«/) 



cosh 27TZ/ — cosh 27T^' 



9-1 

n 

i=i 



l + ^ 2 + ,' 2 ) + j V--' 2 ) 



(72) 
(73) 
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The expression for p§ (xi, X2) can be obtained from (|7lD by the formal substitution 1 for 
the product from s — 1 to s — 2m — 2 and by putting m — 1 in the remaining integral. 

As one of possible applications of the formulae (|7l|)- ([72]) let us consider the dispersion 
of the sizes of localized wave functions. It is seen from ( [71] ) that the distant exponential 
asymptotics of the probability distributions does not fluctuate. On the other hand, it would 
be natural to define the wave packet size i?g as some integral property. For example, let us 
define i?g as: 

4 L r 

R~e 1 = zJ dx \^)\- ( 74 ) 

-L 

The coefficient 4/3 cancels the mean value of the fast oscillating factor sin 4 (fcx + 5) (see 
([0])). Then we have: 

{R ~ El) = ^l dXPE ^ x) = f = W (?5) 
Here p{E) is the density of states at large E: 

P( E) = ±. (76) 

The expectation value of the square of R^ 1 can be found from the correlation function 

(|^(x)^(x')| 4 ): 

L x 

(Re 2 ) = Jdx J dx'p^\x',x). (77) 

Using the explicit expression for p% ,1 \x l ', x): 



L -L 

(2,1), 







and evaluating the integrals over dx' and dv we obtain: 

1 



(78) 



(Re 2 ) « 0.23^. (79) 



With ( [75] ) it gives us the mean square relative dispersion of -R^ 1 : 

(Re 2 ) - (Re 1 ) 2 



(Re) 



0.13. (80) 



5 Mean current in ID mesoscopic ring with the mag- 
netic flux <& 

Let us consider a one dimensional metal ring in transverse magnetic field. The expectation 
value of the current operator for one electron stationary state becomes non-zero and the 
energy receives T -odd term. Then the Fermi levels for left and right directions of the mean 



13 



velocity turn out to be shifted one about another. As a result, a persistent current flows 
along this ring in the ground state |HJ , |L8] . 



For a complete treatement of this physical system one should take into account many- 
body effects such as electron-electron Coulomb interaction (see e.g. [[19]]), fluctuations of the 
chemical potential f2(| etc. I solve here only the one-particle problem and I compute the 



mean current / corresponding to a one-electron state on the Fermi level. There are arguments 
( [pl|,p2|) that the total current is close to / but, of course, further investigations are needed. 

We will assume the ring size 2L to be comparable the mean free path. Then the local- 
ization effects do not lead to the total suppression of /, but / is rendered to be a nontrivial 
function of the magnetic field (see below). (The case of ordered inhomogeneous conductor 
has been considered in the paper [p4fl). 

There exists the gauge by which the wave function of an electron in the ring with the 
magnetic flux $ obeys the boundary condition: 

ip(L) = exp (2tt«$) ip(—L), (81) 

and the Hamiltonian has the previous form ([!]). The mean absolute value of current cor- 
responding to a state with energy E can be represented in the limit (|5|) as follows p3| . 



3n = -2; -w ■ (82) 



{h — c — e — 1, the magnetic flux quantum is equal to 1.) The condition fl8~l|) is nonlocal, 
therefore, the formula (|82"D for / cannot be rewritten in terms of functions like u(x),u(x) 
of Sections 1,2. It has been shown in M, however, that / can be expressed directly via the 



elements of the T-matrix fll4l) . It is worth noting that we have functional representation 
just for them. 

Indeed, by the construction, the matrix T = T(— L,L) satisfies the "unitarity" condi- 
tions: 

<x z TV = T~\detT= 1. (83) 
Therefore we can parametrize in the following way: 

T — { cos hlV as , sumTe^ \ , . 

~ \ sinhTe"^, coshre"^ ) ' ^ ^ 

where T, a s and f3 s are slowly varying real functions of L. The mapping of the initial data 
space at the point x = —L into the space of solutions of the equation (|7|) at the point x = L 
is realized in the basis (u' ± iku) by the transfer-matrix : 



T = exp(ikLa z )T exp(ikLa z ) = 

(85) 

The condition ( |8"TD is equivalent for the matrix to have the eigenvalue e 10 , 9 = 27r$: 

0. 



coshTe i{as+kL) , sinhre i/3s 
sinh Te~ il3s , cosh Te~ i{aa+kL) 



det (t 
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or 



t(E) = coshrcos(a s + kL) = cos9. (86) 

The last equation defines the set of the energy E = k 2 values, and, using the formula (|8~2|) 
for j n we obtain: 



J2$(E-E n )\j n \ =Y,6(E-E n ) 



sin 9 



r'(E) 



5{t{E) -cos 9) \sm9\ 



(87) 



We can get rid of the 5 - function in much the same manner as in the previous section. 
Averaging over the interval AL (1/k <C AL -C I) of the (half-) circumferences L we find: 



k 



I = 2(—5 (t(E) - cos 9)) |sin 9\ 

LJ 



sinh T + sin 9 < 



k 



sin 9 



sin 9 



2 / + °° 

—= ( dp exp ( — /i 2 (sinh 2 T + sin 2 

v 71 " \-L 



(88) 



It is important that sinh 2 T can be rewritten so that the direct analytical continuation from 
the surface £ becomes possible: 



sinh 2 r = (1,0)7VT ^ J 



39) 



Here the sign t denotes, as usually, transposition, and complex conjugation does not enter 
this formula. For sinh 2 T to have the most simple form in terms of the fields p,^ the field 
ip~ should obey zero-valued initial condition: 



<T(-L)=0. 



Then substituting fl2"0p into fl55|) we obtain: 



sinh 2 T = ip~(L) J dt tp + (t) exp ^-i J pdrj 



(90) 



(91) 



The right-hand side of ( pi]) is bilinear in the fields if) but it is non-local. It is convenient to 
start from the transformation of the Hubbard- Stratonovich kind: 



1 = 2 



k sin 9 



f I f If 

/ dp exp (y—p 2 sin 2 9^j / exp < —p 2 ip~(L) / dtip + (t) exp J — i / pdr 



k sin 9 



+oo , 

J dp J dz dz* exp i—p 2 sin 2 9 



\z\ 2 ) X 



x ^exp I — z/j2^> (L) — z/jz* / dtijj + (t) exp ^— i J prfrj jy . 

turning the exponent to the linear in if)^ combination. Retracing the same path 
in the previous section, and introducing the variable £(£) in a manner like ([57]): 



(92) 



as 



£ = -2(1 + 4a)r) + Aiap, 
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Z(L) = 
T>pVr] = const T>pD£ 
we bring the DpDx-integration to the Gaussian form. Performing it we find: 

+00 



(93) 



/ = exp(-^ 



| sin6> 



-jfe N' J d/j f dz dz* exp (—fi 2 sin 9 — \z\ 2 ) x 



J V£exp(-± J dx(e + a 2 ^\z\ 2 e~ 
;)=o V -l v 



Let us change /x for the integration variable cr: 



2 2| |2 

a /i |z| 



(94) 



(95) 



and shift the trajectory by —a: 



-> f (ar) - cr. 

Then the integral over (zz*)-plane is calculated exactly and we obtain the representation of 
I as the matrix element: 



I = exp(— 



k 



sin 9 



(T 2 (0|exp(-2L^)|T 1 (0), 



Here H is defined in (|67|) , and the functions Ti^^) have the following forms: 

T 1 (0 = exp 



(96) 



(97) 



and 



T 2 (0 = exp I-- exp [ -| 



sm( 



(9c 



Using the expansion in terms of the complete set (pq) and the integral representation for 
K2iv{y) we find the explicit formula for the current /: (Z~ x = a/2) 



k (I 



1/2 



— — — exp —— sin 9 



71 L \L 



21 



+00 



dt cosht 



,2/i , „:„u27 ex P I "TTF* 



sin + sinh i 



2L 



(99) 



When sin 2 9 runs from to 1 , the value of I increase monotonically from to I max ■ The 
first order of the Taylor series in 9 agrees with the result of ||. It is worth noting that 
the parameter l/L determines not only the absolute value of /, but its dependence on the 
magnetic field as well: 



sin 2 fl 



+00 



+00 
f v cxp 

J cosh y " 
-00 y 



dt cosht 
sin 2 9 + sinh 2 1 



exp 



I 2 

1 2 

2L 



(100) 



This formula allows us, in principle, to find l/L from the run of the experimental curve 1(9). 
Finishing this section I would like to emphasize that all its formulae are exact in the limit 
(H) and no analogue of the "non-resonance" terms have appeared. 
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6 Conclusion 



Let us suppose for the random potential U(x) to have the finite correlation length k x : 

< U(x)U(x') >= -DKexp{-K\x -x'\) . (101) 
We can take it into account in the limit 1 <C nl by the renormalization of the parameter a: 

^ = TT%5- (102) 

It should be noted that this renormalization may be sufficient since the inequalities 1 C kI 
and ~ 1 can take place simultaneously. 

Indeed, the correlator (|101|) corresponds to the following measure of the integration over 
the fields (±: 




2>C±exp<~ / [l|C| 2 + 2z^(CV-C + C'-) + (l + ^)lCr) dx). (103) 



On performing the "bosonisation" and passing to the variable £ we obtain some effective 
action S e ff(£). The terms with derivatives in the exponent of ( |103| ) would produce the 
terms of S e ff(£) containing derivatives as well. The localization length I is the only parameter 
of the length dimension occurring in the unperturbed problem. Therefore the contributions 
of those " non-mar kovian" terms would be suppressed by powers of the quantity (k/) _1 . 
Neglecting them we come to the formula ( |102| ). 

The variable a' appeared in (|64|) , (|65D could be considered as the global order parameter 
corresponding to the localization. In fact, the non-zero value of the correlator pe{x,x') in 
the thermodynamic limit is the consequence of the following relationship: 

lim — = - > 0. (104) 

On the other hand, the quantity A = exp(— a') is conjugated to the wave function norm. 
Then the inequality (|104 ) corresponds to the exponential in the average increase of the 



functions u, u (|7|). 

It is worth noting that the large-scale behaviour of wave functions is governed by some 
averaged characteristics of the potential U(x). In the perturbation theory framework they 
appear as some infrared singular integrals. The integrands are multipoint products of po- 
tential with fast-oscillating exponentials. Thus, if the quantity 

M..*>=r^(*(,-^ )IW , d05) 



k- 1 < A < / 



x-A 

starting from some A 



becomes independent on x, the Abrikosov-Ryzhkin model with some effective a and a can 
be used to investigate the properties of the wave functions. 

A method similar to that presented here allows one to derive path integral representation 



for any averaged combination of Green functions at arbitrary energy E |25|. Unfortunately. 



I succeeded in computing such path integrals in some of the simplest cases only. 
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